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The effects of heat loss on the structure of laminar flamelets, which are the constitutive elements of turbulent 


flames under the most common operating conditions, are investigated for typical aeronautical gas-turbine operat- 


ing conditions at take-off. The magnitude of heat loss is quantified via the “enthalpy defect” measured with re- 


spect to an adiabatic flame. A procedure to generate laminar flamelets with an assigned enthalpy defect at the 


boundaries is devised and applied to nonpremixed propane/air flames, as propane reproduces the essential fea- 


tures of higher hydrocarbon combustion. It is found, contrary to commonly held beliefs, 


that the enthalpy defect 


has a significant effect on the concentration not only of minor species, but also of main reaction products. Such 


effects are found in general to be more pronounced for fuel-rich conditions. An impact is anticipated on the for- 


mation rate of nitric oxides. The effects of scalar dissipation rate are also discussed. 


Keywords: Turbulent Combustion, Heat Transfer, Laminar Flamelets, Hydrocarbon Combustion, Aeroengines, 


Gas Turbines 


Introduction 


Models for turbulent combustion in industrial applica- 
tions lean heavily on the “laminar flamelet” approach 
[1-5] which, under suitable assumptions which are found 
to hold in most cases, reduces the value of any state 
quantity ø (e.g., density, temperature, concentrations, 
etc.) ina laminar flame to a function of a single quanti- 
ty, the conserved scalar (or mixture fraction) Z, with the 
meaning both of a suitably scaled elemental mass frac- 
tion (or a linear combination of the mass fractions of the 
different elements) and, under the assumption of equal 
diffusivities, of a scaled enthalpy: 

p= (2 (1) 

In turbulent flows, the mean value of the state quantity 
under consideration can be recovered by introducing an 
appropriately defined probability density function (pdf) 
of the conserved scalar, and then convoluting Eq. (1) 
with this pdf. Higher scalar moments, such as variances, 
can also be obtained in a similar way. 

The model can be extended to account for finite-rate 
chemistry (albeit limited to relatively small departures 


from equilibrium) by introducing a second fundamental 
state quantity, the scalar dissipation rate conditioned at 
the flame front (e.g., at stoichiometric, or at the location 
of peak temperature, or peak heat release), denoted as y,. 
In this approach, termed “stretched laminar flamelet” 
[6-9], any state quantity can accordingly be expressed as 
a function of these two quantities: 
p = o (Z, x) (2) 
and evaluation of mean values in turbulent flows re- 
quires the introduction of a joint pdf of Z and y., which is 
still manageable. The functions (2) are evaluated by ded- 
icated codes handling complex chemical kinetics and 
detailed transport for counterflowing jets of fuel and oxi- 
dizer, which allow (once an appropriate definition of the 
conserved scalar is put forth) to recover all state quanti- 
ties as a function of Z. The conditioned scalar dissipation 
rate y. acts as a parameter (specified indirectly by as- 
signing the velocity of the two counterflowing jets), 
which must be varied from near-zero (equilibrium) to the 
quenching limit y,. For y.>y,, the flamelet is extin- 
guished, and state quantities do not depend any longer on 


Xe 


In high-temperature combustion, heat loss effects (by 
radiation in particular) have a direct bearing on state rela- 
tionships; in such circumstances, a third fundamental 
state quantity is required to define state functions in la- 
minar flamelets (incidentally, Z no longer represents en- 
thalpy in this case). The radiant fraction (fraction of 
combustion-generated heat lost by radiation) is some- 
times used [10]; laminar flamelet profiles for the non- 
adiabatic case at hand are then derived by ad hoc adjust- 
ments of the corresponding adiabatic profiles (e.g., see 
[11]). A more rigorous approach [12] involves introduc- 
ing the enthalpy defect (difference between the enthalpy 
of two flamelets, one adiabatic, the other undergoing heat 
transfer, for the same value of Z). Non-adiabatic flame- 
lets can be calculated by enforcing appropriate boundary 
conditions (see next Section), and lead to the form 

9= 02% 9 3) 
which requires a three-variable pdf for evaluating 
mean values in turbulent flames. Adopting the three- 
variable approach accordingly necessitates a very exten- 
sive flamelet “library”, with both y. and ¢ as a parameter. 
We refer to the whole of the flamelets with a given en- 
thalpy defect as to a “shelf”, and in particular, for ¢ = 0, 
as to the “adiabatic shelf’. In order to get an accurate 
solution, each shelf may require from 10 to 20 flamelets 
(with different values of y.), and in highly radiating 
flames, as many as 15 shelves are used [13, 14]. This 
means that, on the whole, computation of quite a few 
hundred flamelets is required. Although the computation 
is performed off-line, this involves a substantial amount 
of computational work. 

This is most often avoided in models of current usage 
by invoking independence of species concentration on 
the enthalpy defect (or, somewhat equivalently, on the 
radiant fraction), and consequently assuming that heat 
transfer affects only temperature (and density). This as- 
sumption is generally thought to be well verified for ma- 
jor species, but in a study [15] concerning a syngas/air 
flame (syngas is a mixture of 40% CO, 30% Hb», 30% N2 
in volume), the authors find that concentration of some 
minor species is significantly affected by the enthalpy 
defect (despite the fact that relatively limited enthalpy 
defects, down to -80 kJ/kg only, are considered there). 

Moving on these grounds, a more ambitious investiga- 
tion is carried out in the present paper, aimed at generat- 
ing a flamelet library to be used to model gas-turbine 
combustion at typical take-off conditions. In this case, air 
is preheated to a high temperature in the compressor, 
while the fuel is liquid kerosene. However, owing to the 
high air temperature, fuel burns effectively very much 
like a vapour [16] (at least at full power operation), then 
phase transition effects can be neglected. Further, higher 
hydrocarbon chemistry is extremely complex, and a un- 


iversally accepted reaction mechanism does not exist. 
Nonetheless, it is well known that propane closely re- 
produces the main features of higher hydrocarbon com- 
bustion, and kerosene in particular [17]. Accordingly, 
fuel is modelled as propane in this study, and typical gas- 
turbine inlet conditions at take-off (Ty; = 852 K, p = 
2.985 MPa) are considered. The main focus of the pre- 
sent paper is on results concerning flamelets subjected to 
low strain (corresponding to a scalar dissipation rate y. ~ 
0.4 s”), but some results are also reported relative to a 
more extensive investigation concerning flamelets with 
strain rates up to extinction, though limited to the adia- 
batic shelf for the sake of conciseness. 

In the next Section, the computational procedure is 
concisely outlined, whereas results are reported in the 
subsequent Section. 


Computational Procedure 


Flamelets are computed by means of the OPPDIF 
code [18], with detailed transport properties and thermal 
diffusion included. The reaction mechanism is taken after 
Peters and Rogg [19], including 34 species (up to C3-) 
and 87 reversible reactions. 

In order to account for heat transfer, the flamelet li- 
brary is organized in “shelves”. The “adiabatic” shelf is 
computed by enforcing temperature boundary conditions 
corresponding to the nominal temperatures of oxidizer 
and fuel (Ta = 852 K, Taer = 300 K), at the respective 
sides. In order to account for heat losses, nonadiabatic 
shelves are introduced, with lower boundary tempera- 
tures with respect to the adiabatic reference shelf. If the 
assumption of unity Lewis numbers for all species is re- 
tained, it is possible to compute a shelf with a uniform 
enthalpy defect by merely assigning appropriate boun- 
dary temperatures. However, when differential diffusion 
is introduced, assigning appropriate temperatures at the 
extremes does no longer allow a uniform enthalpy defect 
across the flamelet. In the following we therefore adopt 
the enthalpy defect at the boundaries (oxidizer and fuel 
side), denoted as ¢;, as the parameter to quantify the ef- 
fect of heat losses, with the understanding however that 
the enthalpy defect ¢ across the flamelet will not be uni- 
form at the value ¢; (see also Fig. 1 below). 

Assume that shelf denoted by index m (i.e., with a 
specified value of ¢) has already been computed, and that 
another shelf m+J is going to be computed, featuring an 
additional enthalpy defect AC (negative) with respect to 
the former. The boundary conditions (temperature and 
species concentration) at the oxidizer and fuel sides must 
be derived accordingly. An immediate way to enforce the 
requested enthalpy defect at the extremes is to lower the 
corresponding temperatures, so that the oxidizer and fuel 


enthalpies exhibit the requested defect ¢+AC with respect 
to the adiabatic shelf. This procedure however is limited 
to small enthalpy defects, because as the temperature is 
reduced below ambient the gases soon condense, or an- 
yway lose their perfect-gas behaviour. An alternative, 
more general procedure is therefore required. 

To this end, we observe that, for a given value of the 
conserved scalar Z, the temperature drops rapidly as the 
enthalpy is reduced, while the species concentrations 
experience only relatively minor changes close to the 
oxidizer and fuel sides (while significant differences may 
arise closer to stoichiometric). With this in mind, the fol- 
lowing procedure is suggested. The mixture enthalpy is 
related to temperature and composition as 

h = SYA; (D) (4) 

If the composition of the mixture of N species is 
known, and the enthalpy is assigned, this can be easily 
inverted to recover the corresponding temperature 

T=T(h, Y%, Yo, ...,5 Yn) (5) 

Points of shelf m close to the boundaries can be 
processed after Eq. (5), with enthalpy assigned as the 
relevant enthalpy of shelf m, diminished by AC, and with 
fixed composition. The first point, moving away from the 
boundary, for which the resulting temperature turns out to 
be equal or greater than a prescribed temperature T“, 
e.g., the ambient temperature (incidentally, different T” s 
might be assigned at the two sides, if required), is chosen 
to enforce the boundary condition for shelf m+/1. 

When detailed transport properties are introduced, 
further problems arise in the definition of the flamelet 
library. While dedicated computational codes compute 
flamelets in physical space, these must be converted in 
conserved scalar space in order to be used by turbulent 
combustion models. Unfortunately, the definition of con- 
served scalar is ambiguous in the presence of differential 
diffusion. The problem is tackled in the literature by in- 
troducing heuristic definitions of the conserved scalar, 
aimed at obtaining a monotonic trend of the conserved 
scalar Z vs. the physical coordinate x. In the present study, 
Z is defined in terms of a scaled mass fraction of element 
nitrogen, Zy: 

Z = (Zn = Zy.air) / (Zx, fuel 7 Zy, air) (6) 

where Zy, air and Zy, juei denote the nitrogen elemental 

mass fractions in the oxidizer and fuel stream, respec- 

tively (the latter being zero in the present case). This is 
found to produce the desired monotonic behavior. 


Results 


The procedure described above to assign the enthalpy 
defect for the shelf under consideration ensures that the 
flamelet features the required enthalpy defect at the 
boundaries, i.e., at the oxidizer and fuel sides. However, 
due to the effect of differential diffusion, the enthalpy 


defect will not be uniform across the flamelet thickness, 
as indicated in Fig. 1. It is in principle possible to com- 
pute the properties of a ‘pseudo-flamelet’ with uniform 
enthalpy defect by interpolating between adjacent shelves; 
this would make the application to turbulent combustion 
models straightforward. The discretization error inherent 
in the interpolation procedure can be reduced to an arbi- 
trarily small proportion by adopting a suitably small AC 
between adiacent shelves. The interpolation procedure is 
important in applications to turbulent flames, but repre- 
sents anyway a sort of higher-order correction to the 
flamelet profiles. Accordingly, it is not considered in the 
following comments, mainly focused at emphasizing that 
there is an effect of the enthalpy level on species concen- 
trations. The range of enthalpy defects at the boundaries, 
¢p, considered in this test case is 0 (adiabatic), -50 kJ/kg, 
=100 kJ/kg, -150 kJ/kg and —200 kJ/kg. 
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Fig. 1 Enthalpy defect, temperature, density: solid line ¢,=0, 
long-dashed line ¢,= -50 kJ/kg, short-dashed line (= 
-100 kJ/kg, dash-dotted line G= -150 kJ/kg, dash- 


double dotted line (= —200 kJ/kg. 


Incidentally, even larger enthalpy defects may take 
place in furnaces, thereby further stressing the need to 
correctly account for the effects under consideration. 

In analyzing the results below, it must be borne in 
mind that stoichiometric conditions for propane/air 
combustion corresponds to the value Z,, = 0.07. 

Figure | also shows the effect of the enthalpy defect 
on temperature profiles, as well as the ensuing effect on 
density (although the latter is also affected by the mixture 
molecular weight, varying with composition). Tempera- 
ture is of course reduced when increasingly negative en- 
thalpy defects are considered, with the effect being more 
apparent at the fuel boundary. For instance, the peak 
temperature of the adiabatic flamelet is 2562 K, which is 
reduced to 2464 for ¢, = -200 kJ/kg. It is also noted that 
the range of Z spanned by the different curves is reduced 


with increasingly negative ¢,s, as a result of the proce- 
dure described in the previous Section to avoid unrealis- 
tically low temperatures at the boundaries. The high nu- 
merical values for the density are due to the high pressure 
conditions under consideration. 

Figure 2 reports the concentration of major chemical 
species, in terms of mass fractions. In particular, the en- 
thalpy defect is seen to have a major impact on molecular 
hydrogen H2. While its peak concentration is nearly un- 
affected, flames with highly negative enthalpy defects 
feature a major increase of H concentration near the fuel 
side, with the appearance of a secondary peak which 
eventually becomes even higher than the one of the 
adiabatic flamelet. 

Molecular oxygen O, concentration is seen to increase 
with increasingly negative enthalpy defects, indicating 
that combustion is less complete. For instance, if we fo- 
cus on O; concentration at Z=0.02, it is found that for ¢, 
= -200 kJ/kg this is increased by 13% relative to the 
adiabatic flamelet. No O2 is found beyond stoichiometric 
(but this is no longer the case for higher ges, see Fig. 4 
below). 

Concentration of propane C3Hg is somewhat affected 
by the enthalpy defect, even though the effect is not ap- 
parent due to the fact that the ordinate scale spans the full 
range zero to unity. Notice that no overlap is seen be- 
tween O, and C3Hg at the low value of ye under consider- 
ation; the picture is different for higher yes, i.e., increas- 
ing finite-rate chemistry effects, in which case some 
overlap is apparent, similar to what is shown in Fig. 4 
below for molecular oxgen. 


An important effect is observed on water gas H3O, 
which is the major reaction product (together with CO2, 
and eventually CO). Peak concentration is seen to in- 
crease about 15% with respect to adiabatic for ¢, = -200 


kJ/kg, thus underlining the need to accurately account for 
enthalpy defect effects in combustion models. 

Another important species featuring a marked effect of 
the enthalpy defect is carbon monoxide CO, with the 
peak value decreasing about 25% (with respect to adia- 
batic) in the range mentioned above. It is anyway seen 
that CO concentration is important only for mixtures op- 
erating near stoichiometric and fuel-rich, which is not the 
case for gas turbines. It is anyway safe to say that CO 
predictions are subjected to considerable uncertainty. 

The resulting effect on CO, turns out instead to be 
much less significant, in the range of about 2% (with 
respect to adiabatic). 

Methane CH, represents (together with other species, 
see below) unburnt hydrocarbons; its concentration is 
seen to be somewhat affected by the enthalpy defect for 
Z>0.17. 

Acetylene C,H, is important both as an unburnt hy- 
drocarbon, and a precursor of soot formation. It appears 
in significant concentrations in fuel-rich combustion, and 
it is seen to show some dependence on the enthalpy de- 
fect for Z>0.13. 

Another unburnt hydrocarbon is ethylene C2H4, which 
is seen to be relatively sensitive to the enthalpy defect for 
Z>0.18. 

Concentrations of some minor chemical species are 
reported in Fig. 3. As they occur in appreciable propor- 
tions only near stoichiometric, the range of the conserved 
scalar reported in the Figures is in this case restricted to 
[0, 0.2]. The enthalpy defect is seen to produce an appre- 
ciable effect on atomic hydrogen H concentration, with a 
peak value reduced by 30% on a relative basis when pas- 
sing from the adiabatic flamelet to the ¢, = -200 kJ/kg 
one. 
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Fig. 2 Concentration of major species: solid line ¢,=0, long-dashed line ¢,= -50 kJ/kg, short-dashed line (= -100 kJ/kg, 
dash-dotted line (= -150kJ/kg, dash-double dotted line (= —200 kJ/kg. 


Fig. 3 Concentration of minor species: solid line ¢,=0, long- 
dashed line = -50 kJ/kg, short-dashed line ¢,= -100 
kJ/kg, dash-dotted line G= -150kJ/kg, dash-double 
dotted line €,= -200 kJ/kg. 


A particularly important effect is noted for atomic 
oxygen O, with peak concentration decreasing by over 
35% in the same range as above. This species is relevant 
to nitric oxide emissions, as the production term due to 
the thermal mechanism can be expressed as 

z a = =| 2 Xo Yva 

Wyoth = 2M yoB,T exp( R P M, My (7) 

where the M;s indicate molar masses, and rate con- 
stants refer to the rate-limiting reaction 

N,+O—>4NO+N (8) 

The effect of reduced oxygen concentration with in- 
creasingly negative ¢ adds up to the accompanying tem- 
perature decrease. It is therefore expected that taking into 
account the effect of enthalpy defect can lead to refined 
predictions of nitrogen oxide emissions from flames and 
combustors. 

Hydroxil OH, which affects soot oxidation, is another 
rate-limited species, and is seen to be significantly af- 
fected by the enthalpy defect (peak concentration reduced 
by up to 20% in the range quoted). Anyway, it is also 
very sensitive to yc. 


The reaction intermediate hydrogen peroxide H,O, 
also shows a significant effect of ¢ (peak reduced by over 
20%). 

Methylidyne CH is closely associated with formation 
of prompt NO; the enthalpy defect appears to bring about 
relatively minor effects on its concentrations. At any rate, 
level are so small that significant formation of prompt 
NO seems unlikely (though this might be an effect of the 
approximation involved in identifying fuel as propane). 

Formyl CHO is seen to occur in trace concentration 
only. Anyway, it is often used as a reaction zone marker 
in hydrocarbon combustion; its concentration is observed 
to be very sensitive to enthalpy defect, as the peak is re- 
duced by over 40% when passing from adiabatic to ¢,= 
-200 kJ/kg. 

Lastly, some effects of the impact of finite-rate chemi- 
stry, as identified by the conditioned scalar dissipation 
rate y,, are reported in Fig. 4 for the adiabatic shelf only. 
As many as 17 flamelets, for stretch conditions ranging 
from near-equilibrium to near-extinction (which is found 
to occur beyond y= 327 s") have been generated, though 
only 5, plus the inert mixing state, are shown in the Fig- 
ure in order to avoid clutter. It is seen that while in 
near-equilibrium conditions the mass fraction of molecu- 
lar oxygen is basically zero beyond stoichiometric, for 
increasing es some oxygen can survive on the rich side 
of the flame due to insufficient residence time for com- 
plete combustion. Further, acetylene concentration is 
shown to be quite sensitive to stretch conditions, with an 
attendant significant effect on soot formation. Similar 
significant effects are found for the concentration of 
atomic oxygen, appearing in Eq. (7) for the nitric oxide 
formation rate, and of hydroxyl, an important soot- 
oxidizing agent. 
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Fig. 4 Concentration of molecular oxygen and acetylene: 
solid line y= 0.4 s"', long-dashed line y. = 82 s“, 
short-dashed line y. = 278 s, dash-dotted line y, = 
309 s, dash-double dotted line No = 327 s! , dotted 


line inert state. 


Conclusions 


Present results emphasize the need to correctly ac- 
count for the effect of heat loss on the concentration of 
chemical species in laminar flamelets, which in turn are 
the thermochemical database passed as an input to turbu- 
lent combustion models. Effects are seen to be quite re- 
markable not only for minor species, but also for some 
major reaction products such as H,O and CO. Effects are 
in general more pronounced at the fuel side, though sig- 
nificant effects are also seen for O, which is only present 
at the lean side (at least for the low value of y, consi- 
dered in this study). Further, it is as well to remark that 
the range of enthalpy defects considered in this study is 
quite limited with respect to the one which can actually 
take place in furnaces, where heat transfer to the wall is 
the goal. 

Future work will involve exploration of the feasibility 
and possible advantage of the approach of the ‘pseudo- 
flamelets’ with uniform enthalpy defect ¢, to be generat- 
ed by interpolation, and the evaluation of the effect of 
heat losses on rate-limited species, such as nitric oxide. 
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